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Abstract. The Atmospheric Chemistry Experiment ACE satellite (SCISAT-1) was 
launched into an inclined orbit on 12 August 2003 and is now recording high signal-to- 
noise 0.02 cm' 1 resolution solar absorption spectra covering 750-4400 cm' 1 (2.3-13 pm). 
A procedure has been developed for retrieving average dry air C0 2 mole fractions (X C 02 ) 
in the altitude range 7-10 km from the SCISAT-1 spectra. Using the N 2 continuum 
absorption in a window region near 2500 cm' 1 , altitude shifts are applied to the tangent 
heights retrieved in version 2.2 SCISAT-1 processing, while cloudy or aerosol-impacted 
measurements are eliminated. Monthly-mean Xco 2 covering 60°S to 60°N latitude for 
February 2004 to March 2008 has been analyzed with consistent trends inferred in both 
hemispheres. The ACE Xco 2 time series have been compared with previously-reported 
surface network measurements, predictions based on upper tropospheric aircraft 
measurements, and space-based measurements. The retrieved Xco 2 from the ACE-FTS 
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spectra are higher on average by a factor of 1.07±0.025 in the northern hemisphere and 
by a factor of 1.09±0.019 on average in the southern hemisphere compared to surface 
station measurements covering the same time span. The ACE derived trend is -0.2% yr' 1 
higher than measured at surface stations during the same observation period. 

Running Title : Carbon Dioxide Trend Measurements 

1. Introduction 

At present there are no global satellite -based measurements of C0 2 profiles despite a 
continuing increase as a function of time with important implications for climate change 
[IPCC, 2007]. Inversion of atmospheric data for the estimation of carbon dioxide (C0 2 ) 
surface sources and sinks relies primarily on the availability of atmospheric remote 
sensing and ground-based surface network observations. However, current observations 
are sparse, particularly over ocean regions, the tropics and subtropics, and in the southern 
hemisphere. Recent studies demonstrate that an average precision of 2.5 ppm (1 
ppm=10' 6 ), less than 0.7% of 375 ppm (the value of C0 2 dry air mole fraction (X C 02 ) 
circa 2009), would be required from monthly averaged column data on a 8 x 10 km 
footprint to achieve comparable constraints on surface sources as compared to those from 
the existing flask network [Rayner et al., 2001], Global column measurements of Xco 2 
with a precision of 1 ppm (0.3%) are required to quantify the variation of sources and 
sinks and improve future climate forecasts [Miller et al., 2007]. Pak and Prather, [2001] 
concluded that satellite observations of upper tropospheric C0 2 profiles have the 
potential to provide a constraint for deriving net carbon fluxes, although a precision of 
about 1% is needed. That work suggested a FTS (Fourier Transfonn Spectrometer) 
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would offer the potential for such measurements, taking advantage of high signal-to-noise 
and also simultaneous profile measurements of a large number of additional species down 
to an altitude of about 5 km. 

The first dedicated spacebome remote sensing missions to measure CO? were planned 
to begin in early 2009 with the launch of the Greenhouse Gases Observing Satellite 
(GOSAT) by JAXA (Japan Aerospace Exploration Agency) and NASA’s Orbiting 
Carbon Observatory (OCO). GOSAT was successfully launched on January 23, 2009 
with initial functional testing and validation scheduled to be completed three months after 
launch by the Japanese National Institute of Environmental Studies and the Ministry of 
the Environment (MOE). GOSAT uses a FTS with high optical throughput, spectral 
resolution, and wide spectral coverage. The instrument measures both CO 2 and CH 4 
from thermal infrared and reflected sunlight observations (0.76 to 15 pm) [ Hamazaki et 
al., 2005]. The OCO mission [ Crisp et al. 2004] had a primary goal of measuring precise 
Xco 2 columns in the near infrared from low Earth orbit and quantifying CO 2 sources and 
sinks globally from the measured columns. The OCO mission was designed to use three 
bore-sighted, long-slit nadir-viewing near-infrared grating spectrometers. However, 
OCO failed to reach orbit after its launch on February 24, 2009. 

High precision Xco 2 measurements are required because of the difficulty in 
disentangling the contributions of different sources and sinks [Tans et al., 1990], the 
presence of clouds and/or aerosols, the effects of CO 2 source signal dilution with height, 
and near uniform Xco 2 background profile distant from major source/sink regions. 
Inversions that provide C0 2 fluxes are sensitive to assumptions of transport [Law et al., 
1996; Fan et al. 1998] with currently limited capability in quantifying regional 
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distributions of sources and sinks from atmospheric measurements due to limited 
sampling [ Fan et al 1998; Chevallier et al. 2005; Tiwari et al. 2006]. Transport out of 
the boundary layer to the free troposphere is not well understood and may be key to 
identifying sink regions [ Strow and Hannon, 2008], Knowledge of the present oceanic 
and terrestrial uptake would have immediate implications for a greenhouse gas emission 
policy [Tans et al, 1996]. 

The potential utility of CO 2 tropospheric solar occultation profile measurements with a 
FTS from inclined-orbit measurements for quantifying CO? source distributions has been 
assessed [P atria et al., 2003] assuming a realistic cloud distribution. Those results 
suggest that such observations would be valuable if sufficient precision is achieved and in 
particular would help constrain source flux uncertainties for the tropical lands with no 
existing observations. The Atmospheric Chemistry Experiment (ACE) [Bernath et al, 
2005] has the potential to obtain high precision measurements directly relevant to carbon 
cycle science objectives, given quantitative procedures for cloud and aerosol detection. 
A previous study investigated the feasibility of retrieving CO 2 profiles in the 5-25 km 
altitude range from the ACE-FTS [ Foucher et al, 2009], working with synthetic spectra 
but also providing estimates for retrievals based on real ACE-FTS data. The importance 
of potential bias in CO 2 retrievals due to errors in the temperature and the nitrogen 
pressure-induced coefficients also was noted. 

Seasonal and latitudinal variations of CO 2 mid-tropospheric zonal mean mixing ratios 
[ Chahine et al, 2008] and a 4-year climatology of lower tropospheric clear sky ocean- 
only CO 2 measurements between ±60° latitude [ Strow and Hannon, 2008] have been 
reported recently from nadir AIRS measurements. Measurements were limited to the 
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northern hemisphere and low latitudes in the southern hemisphere. A climatology of 
lower tropospheric SCIAMACHY/ENVISAT Xco 2 total column measurements obtained 
between 2003 and 2005 primarily in the northern hemisphere over land have been used to 
measure seasonal and latitudinal averages and determine an increase rate for that time 
period [ Buchwitz et al., 2007]. 

2. ACE Measurements 

The ACE satellite was launched on 12 August 2003 into a 74° circular orbit at 650 km 
altitude and has been in operation for science measurements since February 2004, 
providing up to 30 solar occultations per day [Bernath et al., 2005]. The ACE orbit is 
optimized for high latitude coverage in both hemispheres. Measurements are recorded 
with the ACE-FTS, a UV-visible spectrometer, and two visible imagers below about 150 
km altitude. The ACE-FTS records spectra at 0.02 cm' 1 spectral resolution (maximum 
optical path difference of ±25 cm) from 750 to 4400 cm' 1 . Full resolution spectra are 
recorded in 2 s with an altitude spacing determined by the scan time, typically 3-4 km, 
varying from 2 km for long, high beta angle (angle between the plane of the orbit and the 
Sun) occultations and up to 6 km for zero beta angle occultations. As a result of 
refraction, tangent altitude spacing is typically less than 2 km in the troposphere. The 
FTS has a circular field of view (FOV) 1.25 mrad (10' radian) in diameter, which 
corresponds to 3-4 km at the limb in the absence of refraction. Photovoltaic detectors 
with near linear response produce a signal-to-noise ratio of greater than 300:1 throughout 
most of the spectral range. Low Sun solar occultation spectra are divided by exo- 
atmospheric solar spectra from the same occultation to remove non-telluric absorptions 
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yielding measurements between about 85°N and 85°S latitude. Middle and low latitudes 
are recorded only over limited time spans because of the high inclination of the satellite’s 
orbit. All occultations within ±60° latitude with the exception of those inside or on the 
edge of the polar vortex have been included in our analysis covering February 2004 to 
March 2008. 

3. Retrievals 

Pressure-temperature profiles for each occultation have been retrieved as part of the 
procedure for analysis of the ACE-FTS spectra. The assumptions for that analysis, the a 
priori profile for Xco 2 , and its trend vs. time are described in the paper by Boone et al. 
[2005], Versions 1.0 through 2.2 processing are described in that paper, where version 
2.2 is the current set of ACE-FTS results. The model described in that paper assumes a 
linear increase in Xco 2 with time based on the equation adopted previously for analysis 
and retrievals from HALOE (Halogen Occultation Experiment) measurements. 

In the version 2.2 ACE-FTS retrievals, Xco 2 in the middle atmosphere is assumed 
constant, fixed to the value calculated from the model. An additional 5.5 ppm is assumed 
in the troposphere. No provision is made for the variation in Xco 2 with location and 
season. At high altitudes (above ~ 60 km), CO 2 volume mixing ratio (VMR) is fitted, 
employing an empirical function to smooth out unphysical oscillations. Pressure and 
temperature are retrieved down to 12 km. Below 12 km, pressure and temperature are 
fixed to the “reanalysis” data from the Canadian Meteorological Centre (CMC). 

Over 12,500 occultations have been acquired by ACE since launch with current ACE 
(version 2.2) profiles detennined assuming HITRAN 2004 spectroscopic parameters 
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[Rothman et al., 2005]. However, a high fraction of the ACE measurements are 
significantly affected by cloud attenuation at middle and upper tropospheric altitudes. 
We developed a cloud identification procedure using measurements in a micro window to 
eliminate those scenes from the ACE database as well as a procedure to improve the 
precision of Xco 2 retrievals from the remaining ACE spectra based on analysis of a 
window with continuum absorption by the infrared fundamental band of N 2 . We next 
describe the procedure for identifying and excluding ACE measurement scenes 
significantly impacted by cloud extinction. 

A single microwindow (956.7-957.0 cm' 1 ) was chosen for the detection of clouds or 
aerosol extinction [Rinsland et al., 1998], As noted in that work, radiative transfer 
calculations with known molecular lines, absorption cross sections, and a reference set of 
vertical gas concentration profiles predict a transmittance of 0.99 and 0.97 above 150 
mbar (-14 km altitude). We analyzed all spectra in the version 2.2 ACE-FTS database at 
60°N to 60° S latitude applying a criterion for each occultation that reaches a tangent 
height of 10 km or below. We assume that if the transmission in the measured spectrum 
divided by that calculated for the simulated spectrum in the window is less than 0.8, then 
the occultation is cloud impacted and not considered for further analysis. Figure 1 
provides a plot of the fraction of remaining measurements vs. latitude after adopting this 
cloud/aerosol rejection criterion. Measurements penetrate on average approximately 5 
km below the tropopause at all latitudes, a result that is consistent with a previous study 
based on SAGE (Stratospheric Aerosol and Gas Experiment) II 1.02 pm observations 
[ Wang et al. 1998], Profiles of Xco 2 were retrieved from the remaining measurements 
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after applying a vertical offset to the tangent altitudes obtained from ACE version 2.2 
processing. We next briefly describe the procedure we adopted. 

As pointed out previously [Foucher et al., 2009], two main difficulties for retrieving 
CO2 from ACE spectra are: (a) accurate determination of the instrument pointing 
parameters (tangent heights) and (b) sensitivity of retrieved CO2 to pressure/temperature 
profiles. Studies by Rinsland et al. [2004] using ATMOS (Atmospheric Trace MOlecule 
Spectroscopy) solar occultation FTS spectra recorded at 0.01 cm' 1 resolution indicated 
that analysis of continuum pressure-induced absorption by N2 and CO2 in the window 
region 2395-2525 cm' 1 could help achieve the goal of precise determination of 
stratospheric densities from occultation measurements. Thus, we have implemented a 
procedure to determine the altitude shift to apply to version 2.2 ACE tangent heights as 
our means to resolve the difficulty (a) mentioned above. 

Continuum absorption by the pressure-induced fundamental band of molecular 
nitrogen at 4.3 pm was studied previously based on FTS measurements [Rinsland et al., 
1989, 2004]. The absorption is significant over a wide altitude range and increases with 
the square of the density. Collision-induced absorption coefficients for N2 (cm 1 atm' 1 ) 
and their dependence with temperature have been reported for 2125 to 2600 cm" 1 
[ Lafferty et al. 1996, Table 1]. Continuum absorption by CO2 also occurs in the window 
region, but its absorption is negligible in the 2500.6-2500.9 cm' 1 interval selected for the 
current analysis [Rinsland et al., 2004, Fig. 2]. The window contains a weak CH4 line at 
2500.7616 cm" 1 with a ~5% line center absorption depth at 6.6 km tangent altitude. We 
determined the vertical offset to apply to individual version 2.2 spectra based on this 
window by comparing measured and calculated transmittances in the 7-12 km altitude 
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range. No vertical offset was applied to the version 2.2 tangent heights above 12 km 
altitude. In the northern hemisphere, the average values of the applied tangent height 
shifts are in the range of -0.19 km to 0.34 km between 7 km and 10 km. For the southern 
hemisphere, they range from -0.25 km to 0.37 km. 

The ACE spectra obtained after the altitude offset is applied were analyzed with a 
multi-microwindow Levenberg-Marquardt nonlinear least-squares global-fitting 
technique [ Carlotti , 1988]. Concentration profiles were retrieved from simultaneous fits 
to a sequence of spectra recorded in the limb-viewing mode. Microwindow regions in all 
of the spectra from a given occultation were fitted simultaneously over a pre-selected 
altitude range with an assumed pressure -temperature profile to derive Xco 2 and 
interfering species based on an assumed set of molecular spectroscopic parameters and a 
set of a priori profiles. Examples of studies using this approach are presented in the 
papers by Rinsland et al. [1987, 1991]. The Levenberg-Marquardt nonlinear least 
squares multi-microwindow analysis approach is the same as is used to retrieve profiles 
from ACE-FTS spectra [ Boone et al, 2005]. 

Figure 2 illustrates fitting results for the five selected wavenumber intervals for the 
sunrise occultation srll055 at 6.48 km used to derive Xco 2 from the ACE spectra. The 
windows span 2615-2635 cm " 1 and were selected for the retrievals based on those used 
for ACE retrievals [ Boone et al., 2005, Table 2] primarily because the lines are relatively 
insensitive to temperature and they can be used for occultations at altitudes as low as 5 
km with minimal molecular interferences. We fitted CH 4 and CO 2 in the three lower 
wavenumber intervals. Only CO 2 was fitted in the two higher wavenumber windows. 
Note that the lines in these windows belong to a subsidiary isotopologue of CO 2 
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( O “C O). There is an inconsistency between these lines and lines from the main 
isotopologue ( 16 0 12 C 16 0), likely a combination of intensity errors for these weak lines 
and atmospheric fractionation. The COO species is a stable isotopologue with the 

18 i6 

0/ O ratio in atmospheric CO 2 influenced primarily by oxygen exchange with water in 
soils and plants [ Francey and Tans, 1987]. Using measurements between 12 and 20 km, 
Boone et al. [2005] determined a discrepancy of 3.5% between these “COO lines and 
lines from the main isotopologue. Intensities for the “COO lines were thus increased 
by 3.5% whenever they were employed in the ACE-FTS pressure/temperature retrievals. 
The value of this “correction factor” has been adjusted to 4.3% for the next processing 
version of the ACE-FTS (version 3.0). 

ACE-FTS spectra are unapodized, but the instrument suffers from a self-apodization 
of unknown origin [ Boone et al., 2005]. In order to fit lines to the noise level in the 
spectra, it is necessary to introduce an empirical function to account for self-apodization 
effects. The analysis here adopts an empirical model that retains the nominal ACE-FTS 
optical path difference (OPD) but multiplies the ACE-FTS modulation function by a 
triangular function that decreases linearly from 1 .0 at zero OPD to a best- fit value at the 
maximum OPD. The form of the empirical function used in the official ACE-FTS 
processing was different [ Boone et al., 2005], mostly because the official processing 
needed to account for the variation with wavenumber of the self-apodization effects. The 
simple function assumed in this work is sufficient to account for instrument self- 
apodization effects over the narrow wavenumber range employed in this study. 

Table 1 provides a list of the most important sources of random and systematic error 
and provides estimates of uncertainties in the ACE-retrieved 7-10 km X C o 2 - Random 
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errors are much less important than systematic errors because of the high signal-to-noise 
ratio of the ACE-FTS spectra. Random errors due to the finite signal-to-noise ratio of the 
ACE spectra is further reduced by averaging the transmittance over the 2500.6-2500.9 
cm" 1 interval, the window region used to estimate the contribution of the N? continuum to 
the spectrum. Random noise is also reduced by reporting monthly mean Xco 2 over 7-10 
km altitude. Temperature is classified as a random source of error, with validation 
studies indicating that agreement of version 2.2 temperatures with other sensors is 
typically better than 2 K in the stratosphere and upper troposphere [Sica et al., 2008]. 
The magnitudes of 7 sources of potential systematic error are considered individually and 
each is discussed below. 

(1) Retrieval algorithm . The algorithm used in the current study was used in similar 
retrievals of ATMOS spectra [Rinsland et al., 1987, 1988, 1991] and ACE spectra 
[Rinsland et al., 2007, 2009] and on the basis of those results and use in comparison 
with model predictions, we estimate the bias in retrieving molecular profiles as 1%. 

(2) COt line intensities. Line intensities for the CCfi bands in the region were derived 
from 0.01 cm" 1 resolution laboratory spectra [ Malathy Devi et al, 1984]. The 
measurements were recorded at room temperature and low pressure in a 6 m base path 
White cell. The gas sample was a 99.995% minimum purity carbon dioxide sample. 
The isotopic composition of the sample was not measured. The intensities from that 
work [ Malathy Devi et al., 1984] are a factor of 1.15 higher than those from an earlier 
lower resolution laboratory study [Hoke and Shaw, 1982, Table 1], though intensities 
from an earlier grating spectrometer band intensities were reported as in agreement 
[Malathy Devi et al., 1984], The line intensities reported earlier [Hoke and Shaw, 1982] 
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were calculated assuming the isotopic abundances in naturally occurring CO 2 , which 
was used for collecting the measurement sample. The above-cited studies suggest 
errors in the assumed intensities may introduce errors of up to 15%. To detennine 
which set of CO 2 line intensities is more accurate, a study was undertaken with solar 
absorption data measured by the Kitt Peak FTS (3 1 .9°N latitude, 1 1 1 .6°W longitude, 
2.09 km altitude). Daily average 2.09-10 km altitude Xco 2 from 1977 to 2008 were 
derived from analysis of the spectral region near 2625 cm' 1 with the SFIT2 version 3.93 
algorithm, which is based on the optimal estimation method [ Rodgers , 1990] modified 
with a semi-empirical implementation [Parris et al., 1992; Connor et al., 1995]. The 
spectroscopic parameters from the HITRAN 2004 database [Rothman et al., 2005] were 
assumed in the retrievals. The comparison shows agreement between the Kitt Peak 
daily mean Xco 2 and monthly-mean measurements from two surface stations (BME at 
32.37°N latitude, 64.65° longitude and WIS at 31.13°N latitude, 34.88° longitude) at 
similar latitude. The ratio of the Kitt Peak monthly means for the 2.09-10 km altitude 
range to the monthly mean measurements from the BME and WIS stations for nearly 
the same time span during December 1995 to June 2005 are 1.019±0.024 and 
1.017±0.02, respectively. A total of 58 monthly means from Kitt Peak based on 
measurements from 5 microwindow with spectral intervals similar to those used for the 
ACE retrievals were compared with 115 monthly means from the BME station and 115 
monthly means from the WIS station. Xco 2 from three windows spanning 3203-3326 
cm' 1 are on average 2.2% lower than those from the 5 windows used for the retrievals 
reported in this work. Based on these comparisons and the offset applied in the version 
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2.2 ACE-FTS pressure/temperature retrievals for these lines, we estimated the 
systematic error in the assumed line intensities as 4%. 

(3) Nitrogen continuum absorption coefficients . Our analysis assumes temperature- 
dependent N 2 collision-induced absorption coefficients [Lafferty et al., 1996] covering 
2125-2600 cm' 1 predicted from measurements recorded in the 0-10 atm and 230-300 K 
pressure ranges, respectively. FTS measurements of pure N 2 recorded at 0.5 cm" 1 
spectral resolution and N 2 absorption coefficient measurements were combined with 
grating spectrometer measurements of N 2 and 0 2 mixtures at 193 to 297 [Menoux et al., 
1993] as both were found to be consistent within their experimental uncertainties. The 
empirical model used to analyze the combined set yielded the predicted absorption 
coefficient (cm 1 atm' 1 ) and its dependence with temperature in increments of 5 cm' 1 . 
The region from 2460 to 2560 cm' 1 is important as it contains window intervals 
dominated by N 2 continuum absorption that provide transmissions potentially useful for 
retrievals of vertical offsets at upper tropospheric and lower stratospheric altitudes. We 
tested the consistency of the altitude shifts deduced from the ACE spectra by comparing 
shifts determined from 10 clear sky occultations covering the altitude range, applying 
the altitude shifts as described previously. Altitude shifts from two microwindows 
(2500.6-2500.9 and 2522.9-2522.95 cm" 1 ) agreed to 0.4 km, but both intervals are in 
regions that are sensitive to errors in the predicted temperature dependence of the N 2 
absorption coefficient [Lafferty et al., 1996, Table 1]. We adopted 3% as the relative 
uncertainty based on these comparisons, the sensitivity due to differences in the 
reported measurement sets [Lafferty et al., 1996, Figs. 5 and 6], possible bias due to 
error in the empirical model adopted to calculate the temperature-dependent absorption 



14 


coefficients [Lafferty et al., 1996, Eqs. 4, 8], and uncertainties in ACE upper 
tropospheric temperatures. More recent FTS integrated band intensity measurements 
[Baranov et al., 2005] derived at temperatures up to 360 K are about 3% lower than 
previous results but about 1% larger than reported by Menoux et al. [1993]. New low 
temperature measurements of N 2 collision-induced absorption coefficients at 5 pm are 
planned (Y. Baranov, National Institute of Standards and Technology, private 
communication, 2009). The tangent height retrieval error (in km) for N 2 continuum 
absorption in the 5-10 km altitude range as a function of microwindow width (cm" 1 ) 
derived assuming the absorption coefficients of Lafferty et al. [1996] yielded an 
uncertainty of 0.035 cm 1 [Foucher et al., 2009, Table 6]. As our primary N 2 window is 
0.3 cm" 1 wide, a retrieval error of 0.3% was estimated. The corresponding uncertainty 
in the Xco 2 is 0.1%. However, this estimate does not take into account the much larger 
spectroscopic uncertainties. We adopted 2% as the estimated uncertainty inXco 2 due to 
uncertainty in the N 2 pressure-induced fundamental band absorption coefficients near 
2500 cm' 1 . 

(4) CO? isotopologue correction. The HITRAN database incorporates assumed 
terrestrial isotopologue abundances for each molecular species a.v defined by De Bierve 
et al. [1984]. We have assumed those values. As the relative abundances of different 
C0 2 isotopologues vary with ecosystem and such differences are not yet well 
understood, but they are recognized as potential indicators for identifying different 
emission sources. Deviations from the assumed HITRAN isotopologue value are a 
potential source of systematic error in the ACE CO? retrievals. As mentioned earlier, 
lines of 16 0 18 0 16 0 are used in the low altitude retrieval of carbon dioxide, and when 
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they are combined with those for 16 0 12 C 18 0 and the 16 0 12 C 16 0 at 12-20 km a bias 
occurs. Based on those results and other studies [ Francey and Tans, 1987; Allison and 
Francey, 2007; Nakazawa eta/., 1997], we estimate atmospheric fractionation among 
different CO 2 stable isotopologues may introduce a systematic error of <1%. 

(5) Aerosols. Extinction due to aerosols contributes to the absorption along the long 
paths measured during solar occultation observations, and in the past significant 
infrared extinction has been measured shortly after major volcanic eruptions (e.g., 
[Rinsland et al ., 1984, Kent et al., 1995]). Currently, aerosol levels are very low, 
and absorption by aerosols is much less significant in the infrared than the visible. 
Additionally, the spectral region for measurement of pressure-induced absorption by 
nitrogen is located in a region of low infrared extinction [Rinsland et al., 1984]. We 
assume an upper limit of 1% as the potential bias introduced by not including 
aerosol extinction in our retrievals. More recently, active remote sensing 
measurements of aerosols and clouds have become available with lidar techniques 
with global distributions [ Winker et al., 2007]. 

(6) Uncertainty in the ACE-FTS ILS. As mentioned earlier the approach used for fitting 

the ACE-FTS ILS differs from that adopted in the standard processing of atmospheric 
spectra, though there is agreement that the measured ILS is broader than predicted 
with the nominal parameters for the ACE ILS. As both studies reach a similar 
conclusion, we estimate an upper limit for the impact of the error in retrieving Xco 2 
at 7-10 km as less than 1%. 

(7) Uncertainty in the altitude shift . Potential bias in detennining the altitude shift is 
linked to systematic errors in the assumed N 2 continuum absorption coefficients 
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(error source #3). The window for detennining the vertical altitude shift is 
relatively narrow (2500.6-2500.9 cm’ 1 ), and we have assumed the altitude shift 
above an altitude of 12 km is zero. Considering only these two factors, we estimate 
a potential error in Xco 2 from both as less than 1%. 

4. Results 

The 7-10 km altitude range for the retrievals was chosen to reach as low as possible but 
retain sufficient observations to allow measurement of statistically-significant seasonal 
changes in Xco 2 and its trend in both hemispheres. The time series of monthly averaged 
measurements was fitted with the expression 

C A = a 0 + ai (t-t 0 ) + a 2 cos (2 n (t-t 0 -cp», ( 1) 

where Ca is the Xco 2 at time t, ao is the mean Xco 2 at time to, ai is the trend, a 2 is the 
amplitude of the seasonal cycle, and (p is the phase corresponding to the seasonal cycle 
maximum. Measurements of Xco 2 and its trend are limited primarily to those from surface 
station sampling sites. Equation 1 is an approximation. 

Figure 3 compares monthly mean ACE Xco 2 at 7-10 km altitude with no altitude shifts 
and monthly average Earth System Research Laboratory (ESRL) NOAA (National 
Oceanic and Atmospheric Administration) surface measurements for the same time 
period. Measurements for 0° to 60°N latitude are in the upper panel and those for 0 to 
60° S latitude are in the lower panels with residuals (measured minus calculated monthly 
averaged Xco 2 ) on an expanded vertical scale. Each individual occultation that results in 
a successful retrieval was used to perform linear interpolation onto 1-km altitude spacing 
before calculating average Xco 2 for the 7-10 km layer. Open triangles in Figure 3 display 
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ACE monthly mean Xco 2 measurements and the vertical lines indicate their standard 
deviations. The regression fit with equation 1 is shown as a solid curve and the slightly 
tilted straight line represents the constant term plus the trend term. Listings of the 
locations and altitudes of the northern and southern hemisphere NOAA/ESRL surface 
station used for comparison with ACE measurements are provided in Tables 2 and 3, 
respectively. The station codes are shown in the figure. 

Figure 4 shows the comparison of the northern (0-60°N latitude) and southern 
hemisphere (0-60°S latitude) monthly mean Xco 2 time series at 7-10 km altitude when 
ACE-FTS tangent heights were adjusted through the analysis with the N? continuum prior 
to the retrievals. As in Figure 3, open triangles represent the measurements and vertical 
lines indicate their 1 -sigma statistical uncertainty with a solid curve indicating the best-fit 
to the time series assuming Equation 1 and residuals displayed above on an expanded 
vertical scale. The time series are again compared with the NOAA/ESRL monthly mean 
time series for the same time period (same as in Figure 3). As can be seen from this 
figure, a significant seasonal variation occurs whereas none is significant in Figure 3, 
reflecting the use of a non-seasonally varying a priori profile below middle atmospheric 
altitudes in the version 2.2 ACE temperature/pressure retrievals [Boone et al., 2005]. 

5. Discussion 

Table 4 provides best fit results to ACE monthly average time series with equation 1 for 
northern hemisphere (0-60°N latitude), the southern hemisphere (0-60° S latitude) and their 
1 -sigma statistical uncertainties. The average X C o 2 , its trend, and corresponding 1 -sigma 
uncertainty obtained from a fit with equation 1 to the NOAA/ESRL monthly mean in situ 
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Xco 2 measurements covering the same time span are also reported. The ACE monthly 
measurements are on average a factor of 1.07±0.025 higher in the northern hemisphere and 
a factor of 1.09±0.019 higher than the NOAA ESRL southern hemisphere surface station 
measurements. 

The mean and standard deviation of ACE Xco 2 for the northern hemisphere covering 44 
months is 412.35±7.20 ppm. The trend and its uncertainty is 3.11±0.57 ppm/yr, 1 sigma. 
The corresponding monthly mean northern hemisphere surface measurement time series 
shows a monthly mean X C o 2 of 382.95±5.98 ppm with a seasonal amplitude of 5.16±0.22 
ppm. The fit to the surface monthly-mean time series from 17 northern hemisphere stations 
indicates a trend of 2.22±0.14 ppm/yr, 1 sigma. The ACE southern hemisphere 7-10 km 
mean and uncertainty of 414.90±6.40 ppm is consistent with the ACE northern hemisphere 
result. The seasonal cycle amplitude and uncertainty given by ACE for southern 
hemisphere is 4.12±1.03 ppm. 

AIRS measurements from very clear ocean-only scenes in both hemispheres between 
±60° latitude over four years with maximum sensitivity at 2-9 km altitude (-300-800 hPa) 
[Strow and Hannon, 2008] show Xco 2 3-9 ppm lower than measured in the boundary layer, 
and therefore the vertically-averaged 7-10 km average Xco 2 likely also should exhibit the 
COt rectifier effect. Analysis shows ACE sampling of minimally-impacted cloud 
occultations and the latitude coverage is too sparse to derive more than one layer with 0°- 
60°N and 0°-60°S latitude averages required to retain a sufficient number of occultations to 
accurately quantify trends and the difference in the seasonal cycles in the two hemispheres. 

We have used the regression fit obtained from the third tenn on the right side of 
Equation 1 to detennine the amplitude of the seasonal cycle for each northern hemisphere 
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surface station and the month corresponding to the seasonal maximum. We limit our 
comparison to the northern hemisphere as the amplitude of the seasonal cycle is small in 
the southern hemisphere [ Ericksen et al., 1996, Figure 3]. The results for the northern 
hemisphere stations are given in Table 5. The month of the seasonal maximum from the 
northern hemisphere surface stations occur between January and April. The amplitude of 
the monthly average seasonal cycle is observed to generally increase with increasing 
latitude (~3 ppm in the subtropics to ~8 ppm at high latitudes) except for the HUN station 
at 46.95°N (amplitude of 13.28±0.56 ppm). The majority of those amplitudes are in the 
range 3 to 6 ppm, consistent with the peak-to-peak seasonal variations of 6-8 ppm 
determined from upper tropospheric northern hemisphere aircraft measurements 
[Nakazawa et al, 1991; Saw a et al., 2008] and ground-based total column measurements 
from Kitt Peak, which showed an average seasonal variation peak-to-peak amplitude of 7 
ppm [ Yang et al, 2002], The 4-year AIRS climatology from ocean-only clear sky scenes 
between 10°N to 60°N latitude with maximum sensitivity to CCT variability at 300-800 
hPa (3-9 km) [Strow and Hanlon, 2008] also reported an average seasonal amplitude of 3 
ppm. 

Figure 5 compares ACE 40°N-60°N latitude monthly average results with results from 
GLOBAL VIEW model weekly time series for 7 ground station locations between 40°N- 
50°N. The fonnat for the ACE measurements is the same as in Figures 3 and 4, and both 
ACE measurements and model predictions cover the same time period (2004 to 2008). 
Equation 1 was used to fit time series of both datasets. The station code for each 
GLOBAL VIEW site [www.esrl.noaa.gov/gmd/ccgg/globalview/] is also displayed. The 
ACE-derived seasonal cycle maximum occurs in the middle of January with an amplitude 
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1.158±0.31%. The month of seasonal cycle maximum for GLOBAL VIEW is found to 
occur in April and the amplitude is 0.92±0.01%. A total of 36 monthly means from ACE 
observations of 217 occultations are included in the analysis. 

The GLOBAL VIEW predictions are derived from NOAA middle and upper 
tropospheric aircraft measurements and are described as having been smoothed, 
interpolated, and extrapolated. The ACE Xco 2 measurements are higher than those 
predicted by GLOBALVIEW, and the phase and amplitude of the seasonal cycle maxima 
as well as the Xco 2 trend are all different. As the GLOBALVIEW database is widely 
used and has been included in validation studies (see for example, [Strow and Hannon, 
2007, Figure 9]), it is an important database, but it is not possible to determine if the 
differences between ACE results and GLOBALVIEW predictions for trend or seasonal 
cycle amplitude and phase are real or artifacts. 

Table 6 provides a listing of ACE monthly average measurements (and standard 
deviations}, a best-fit obtained to the time series with Equation 1, the latitude range, and 
the number of occultations in northern hemisphere for each month from the February 
2004 to February 2008 time period. Table 7 provides a listing of the same information in 
the same fonnat from the southern hemisphere 0°-60°S measurements. It can be noted 
that the fit to the seasonal variation is not as good for the southern hemisphere results as 
those from the northern hemisphere. 

The Xco 2 rectifier effect noted earlier [Strow and Hannon, 2008] indicates Xco 2 in the 
lower mean troposphere is several parts per million lower than in the boundary layer near 
the same site for the same time period. Higher CCL ACE measurements at 7-10 km are 
also inconsistent with model profile calculations for the Mauna Loa station [Ericks en et al., 
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1996, Fig. 1 1] and upper tropospheric aircraft measurements [ Nakazawa et al., 1991; Saw a 
et al., 2008]. 

6, Summary and Conclusions 

Retrievals of Xco 2 have been performed for February 2004 to March 2008 at 0-60°N 
and 0-60°S latitudes from ACE SCISAT-1 spectra based on an analysis procedure that is 
described in Section 3. Monthly average X C o 2 retrieved from the ACE-FTS spectra at 7- 
10 km are higher on average by a factor of 1.07±0.025 in the northern hemisphere and by 
a factor of 1.09±0.019 in the southern hemisphere than NOAA ESRL surface station 
measurements covering the same time span. The ACE derived trend is also ~1 ppm yr' 1 
(0.2% yr' 1 ) higher than measured at those surface stations during the same observation 
period. A bias between remote sensing and correlative measurements has been reported 
previously [ Yang et al., 2002; Bosch et al., 2006], with offsets of several percent between 
the different measurement sets identified. Different biases were reported for the two C0 2 
channels used for retrievals by the AIRS instrument [ Strow and Hannon, 2008]. 
Therefore, there is a need for further laboratory studies of CO 2 in both the visible and 
infrared to reduce systematic sources of error. The systematic errors estimated for our 
results are much larger than those required to quantify the variation of sources and sinks 
and improve future climate forecasts. The results are also not accurate enough to infer 
surface fluxes given the high precision needed for space-based measurements [Miller et 
al., 2007]. Additional low temperature measurements are planned and may reduce the 
uncertainty in the N 2 continuum absorption coefficients, including the atmospheric 
window region near 2500 cm’ 1 . 
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The ACE Xco 2 measurements is inconsistent with the CO 2 rectifier effect reported 
from a 4-year climatology of AIRS lower tropospheric clear sky ocean-only satellite 
measurements [ Strow and Hannon, 2008] and a decrease with altitude in Xco 2 predicted 
above Mauna Loa [. Ericks en et al., 1996]. ACE measurements at 40°N-60°N latitude 
have been compared with NOAA-ESRL surface station measurements for a 36 month 
time period and exhibit agreement in the time of the seasonal cycle peak within one 
month. ACE Xco 2 measurements have also been compared with values predicted for the 
same time span based on the GLOBAL VIEW aircraft database. The phase and amplitude 
of the seasonal cycle as well as the trend for Xco 2 differ from those measured by ACE. 
As indicated in Figure 5, the seasonal amplitudes from ACE and GLOBAL VIEW time 
series are found to be 4.78 ppm and 3.51 ppm, respectively. The time of occurrence of 
seasonal maximum is January and April, respectively. Since GLOBAL VIEW predictions 
represent temporal variations of localized regions (mostly over the U.S.A) and the ACE 
results represent an average over the latitude range 40°N to 60°N, certain differences are 
anticipated. The averaged Xco 2 amount from ACE is higher than GLOBAL VIEW by a 
factor of 1.08. 

Our use of monthly average measurements averaged over latitude range and a 3 -km 
wide altitude range and the need to exclude occultations impacted by clouds, aerosols 
and/or descent reflect the measurement limitations from an occultation instrument such as 
ACE. Despite the sparse sampling, the capability of the ACE-FTS to also precisely and 
simultaneously measure a comprehensive set of chemical tracers with a range of lifetimes 
and sources such as CO, CH 4 , HCOOH, 0 3 , C 2 H 6 , C 2 H 4 , HCN, CH 3 OH, SF 6 , OCS, and 
C 2 H 2 offers the potential to use those measurements over land and ocean scenes in both 
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hemispheres to differentiate air mass types and offer constraints on the factors needed in 
model vertical transport from surface source regions to the upper troposphere [Law et al . , 
1996] (e.g., biomass burning, fossil fuel emissions, continental outflow as a function of 
latitude, season, and time). The ACE near global mid-tropospheric measurements of 
Xco 2 and CH 4 will provide validation for those obtained by the GOSAT FTS [ Hamazaki 
et al., 2005]. ACE observations are just one of many datasets that are needed to help 
quantify the global atmospheric distribution of X C o 2 , to constrain uncertainties in its 
atmospheric flux budget, and improve modeling of lower to upper tropospheric transport. 
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Table 1 . Random and Systematic Sources of Error in Retrieval of 7- 1 0 km Xco 2 


Error Source 

Error 

Type 

Relative Uncertainty (%) 

Finite signal-to-noise 

R 

<1 

Temperature Profile 

R 

1 

Retrieval Algorithm 

S 

1 

Nitrogen continuum absorption coefficients 

S 

2 

CO 2 Line Intensities 

s 

4 

Aerosol 

s 

<1 

CO 2 Isotopologue correction 

s 

<1 

Uncertainty in the ACE-FTS ILS 

s 

<1 

Uncertainty in the vertical altitude shift 

s 

<1 


*S=systematic, R= random, ILS is the instrument line shape function. 
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Table 2. Station codes, locations, and altitudes of the northern hemisphere NOAA/ESRL 
surface stations used in the study 


Station 

Latitude (°) 

Longitude (°) 

Altitude (m) 

BAL 

55.50°N 

16.67° 

7.0 

BME 

32.37°N 

-64.65° 

30 

BSC 

44.17°N 

28.68° 

3 

GMI 

13.43° N 

144.78° 

2 

HUN 

46.95°N 

16.65° 

344 

IZO 

28.30°N 

-16.48° 

2300 

KZD 

44.45°N 

75.57° 

412 

KZM 

43.25°N 

77.88° 

2519 

MLO 

19.53°N 

-155.58° 

3397 

PTA 

38.95°N 

-123.73° 

55 

RPB 

13.17°N 

-59.45° 

3 

SHM 

52.72° N 

174.10° 

40 

THD 

41.05°N 

-124.15° 

107 

UTA 

39.90°N 

-113.72° 

1320 

UUM 

44.45°N 

111.10° 

914 

WIS 

31.13°N 

34.88° 

400 

WLG 

36.27°N 

100.92° 

3810 
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Table 3. Station codes locations, and altitudes of the southern hemisphere NOAA/ESRL 
surface stations used in the study 


Station 

Latitude (°) 

Longitude (°) 

Altitude (m) 

ASC 

7.92°S 

-14.42° 

54 

BHD 

41.42°S 

174.87° 

80 

BKT 

0.20°S 

100.32° 

865 

CGO 

40.68°S 

144.68° 

97 

CRZ 

46.45° S 

51.85° 

120 

EIC 

27.15° S 

-109.45° 

50 

SEY 

4.67° S 

55.17° 

3.0 

SMO 

14.25° S 

-170.57° 

42 

TDF 

54.87° S 

-68.48° 

20 
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Table 4. Best fit parameters to monthly mean ACE northern hemisphere (0-60°N) 
latitude 7-10 km altitude (A), ACE southern hemisphere (0-60°S) 7-10 km altitude (B), 
NOAA/ESRL northern hemisphere surface (C), and southern hemisphere surface (D) 
CO 2 time series for the February 2004-March 2008 time period 


Case 

Mean VMR 
(ppm) 

VMR 

trend 

(ppm/yr) 

Trend 
(% yr' 1 ) 

Seasonal 

Cycle 

Amplitude 

(ppm) 

NMO 

NMS 

A 

412.35E7.20 

3. 1 1+0.56 

0.77E0.14 

5.60±0.92 

44 

333 

B 

414.90E6.40 

2.75E0.61 

0.67± 0.15 

4. 12+1 .03 

46 

455 

C 

382.95E5.98 

2.22E0.14 

0.59±0.04 

5.16±0.22 



D 

378.05±2.62 

1.87E0.07 

0.50±0.02 

0.10±0.09 




Notes: Altitude-shifted monthly mean ACE occultations fit with Equation 1 . NMO and 
NMS are the number of measurement months and the number of valid occultations after 
excluding cloudy and/or aerosol impacted ACE measurements, respectively. Monthly 
mean NOAA/ESRL northern hemisphere (17 stations) (C) and southern hemisphere (9 
stations) (D) were used to infer surface trends and seasonal cycle amplitudes. 
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Table 5. Best- fit for the seasonal cycle amplitude and the month of the peak derived 


from northern hemisphere surface stations 


Station Latitude 

Seasonal Cycle 
Amplitude (ppm) 

Month of Peak 
For Seasonal 
Cycle 

BAI (55.5°N) 

7.82 ±0.61 

January 

SHM (52.72°N) 

7.47 ±0.58 

February 

HUN (46.95°N) 

13.28 ±0.56 

January 

KZD (44.45°N) 

7.69 ± 0.44 

January 

UUM (44.45°N) 

6.71 ±0.50 

February 

BSC (44.17°N) 

6.08 ±0.85 

January 

KZM (43.25°N) 

6.21 ±0.37 

January 

THD (41.05°N) 

3.95 ±0.63 

February 

UTA (39.90°N) 

4.38 ±0.30 

February 

PTA (38.95°N) 

4.82 ±0.65 

February 

WLG (36.27°N) 

4.38 ±0.33 

February 

BME (32.37°N) 

5.55 ±0.31 

February 

WIS (31.13°N) 

4.61 ±0.28 

February 

IZO (28.30°N) 

3.47 ±0.20 

March 

MLO (19.53°N) 

3.29 ±0.14 

April 

GMI (13.43°N) 

3.06 ±0.13 

April 

RPB (13.17°N) 

3.63 ±0.16 

March 


Notes: Date corresponding to the seasonal maximum was determined from a fit based on 
Equation 1. 
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Table 6 . Measured, best-fit calculated, latitude range, number of occultations from ACE 
northern hemisphere monthly averaged measurements between February 2004 and 
February 2008 


Date 

(year month) 

C0 2 (ppm) 
Mon. Avg ± 
Std. Dev. 

Regression 
Fit C0 2 

Latitude range 

NMOC 

2004 Feb. 

409.25 ± 9.79 

410.58 

8.3 - 58.78°N 

11 

2004 Mar. 

406.56 ±6.42 

407.35 

35.92 -59.46°N 

23 

2004 Apr. 

402.78 ±3.98 

406.85 

19.75 -35.64°N 

9 

2004 May 

402.35 ±0.0 

402.42 

57.95°N 

1 

2004 Jun. 

407.77 ± 6.29 

401.55 

51.16 - 51.72°N 

5 

2004 Jul. 

401.10 ± 5.57 

401.78 

33.9 - 58.33°N 

25 

2004 Aug. 

398.48 ±5.42 

403.21 

4.34 - 43.35°N 

19 

2004 Sep. 

404.57 ±4.99 

404.23 

44.47 _ 56.34°N 

11 

2004 Nov. 

414.15 ±8.71 

412.49 

53.88 -59.26°N 

7 

2004 Dec. 

418.26 ±5.40 

413.77 

41.37 -55.44°N 

10 

2005 Feb. 

402.62 ± 
11.44 

413.34 

8.89 - 58.71°N 

3 

2005 Mar. 

409.28 ±7.88 

410.36 

42.42 - 58.24°N 

12 

2005 Apr. 

404.82 ±5.85 

408.80 

4.87 - 32.92°N 

17 

2005 May 

412.76 ±4.57 

406.84 

35.25 -59.25°N 

9 

2005 Jun. 

412.48 ±0.0 

404.57 

44.08°N 

1 

2005 Sep. 

402.13 ±0.37 

407.36 

47.07 - 47.69°N 

2 

2005 Oct. 

413.47 ±9.85 

412.29 

17.88 -36.35°N 

5 

2005 Nov. 

416.54 ±5.25 

415.60 

55.53 -59.94°N 

7 

2005 Dec. 

420.07 ±4.49 

415.92 

50.42 -54.53°N 

3 

2006 Jan. 

419.86 ±7.38 

417.27 

56.57 -59.28°N 

7 

2006 Feb. 

412.70 ± 
11.58 

416.68 

4.79 - 58.69°N 

7 

2006 Mar. 

415.80 ±3.59 

413.61 

54.95 - 59.96°N 

4 

2006 Apr. 

409.90 ±2.31 

412.05 

17.31 -39.38°N 

11 

2006 May 

415.66 ±2.37 

410.06 

30.68 -58.98°N 

5 

2006 Jun. 

409.98 ±0.0 

408.46 

54.76°N 

1 

2006 Jul. 

410.50 ±5.01 

407.72 

48.57 -55.44°N 

5 

2006 Aug. 

406.99 ±7.08 

410.09 

9.98 - 35.54°N 

2 

2006 Sep. 

411. 14 ± 1.91 

410.72 

43.15 -56.67°N 

7 

2006 Oct. 

407.03 ±0.0 

414.53 

24.03°N 

1 
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2006 Nov. 

424.37 ±3.98 

418.76 

56.71 -59.25°N 

5 

2006 Dec. 

423.25 ±3.30 

419.01 

53.32 -56.49°N 

6 

2007 Jan. 

423.99 ±3.97 

420.39 

53.91 -59.91°N 

9 

2007 Feb. 

414.05 ±5.26 

419.72 

1.47 - 57.84°N 

12 

2007 Apr. 

415.39 ± 6.18 

415.39 

8.06 — 41.25°N 

7 

2007 May 

414.14 ± 1.41 

413.26 

26.94 - 58.90°N 

6 

2007 Jun. 

411.96 ±6.36 

411.00 

38.93 -57.64 

4 

2007 Jul. 

410.09 ±0.94 

410.87 

49.37 - 52.23°N 

2 

2007 Aug. 

408.24 ±4.36 

411.715 

17.52 -37.83°N 

7 

2007 Sep. 

410.33 ±2.73 

413.65 

34.35 -55.05°N 

11 

2007 Oct. 

41 1.83 ± 0.0 

418.67 

26.15°N 

1 

2007 Nov. 

427.04 ± 1.38 

421.85 

58.46 -59.97°N 

3 

2007 Dec. 

423.83 ±5.66 

422.69 

47.35 -56.99°N 

10 

2008 Jan. 

427.75 ±6.68 

423.50 

52.55 - 59.78°N 

18 

2008 Feb. 

418.32 ±7.77 

422.63 

22.45 - 51.12°N 

2 


Notes: The best-fit calculated CCfi VMRs were obtained from a fit to the time series 
with Equation 1 . NMOC is the number of number of valid occultations after excluding 
cloudy and/or aerosol impacted ACE measurements and those inside or on the edge of the 
polar vortex. See text for details. 
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Table 7. Measured, best-fit calculated, latitude range, number of occultations from 
ACE southern hemisphere monthly averaged measurements between February 2004 and 
March 2008 


Date (year 
month) 

C0 2 (ppm) 
Mon. Avg ± 
Std. Dev. (ppm) 

Regression Fit 
C0 2 

Latitude range 

NOCC 

2004 Feb. 

401.34 ±0.0 

405.323 

23.44°S 

1 

2004 May 

422.05 ±0.91 

412.867 

58.19 -58.94°S 

4 

2004 Jul. 

419.73 ±3.79 

414.490 

55.80 -59.22°S 

6 

2004 Aug. 

407.94 ±7.81 

414.431 

1.44 - 21.48°S 

4 

2004 Sep. 

416.71 ±3.84 

412.447 

43.04 -59.99°S 

14 

2004 Oct. 

406.25 ±3.57 

411.095 

9.31 -38.74°S 

24 

2004 Nov. 

413.67 ±4.70 

409.847 

39.24 -59.91°S 

42 

2004 Dec. 

416.07 ±8.66 

408.274 

43.64 -52.51°S 

6 

2005 Jan. 

406.54 ±5.53 

407.946 

27.92- 58. 12°S 

22 

2005 Feb. 

405.80 ±4.25 

408.882 

14.74 -40.03°S 

10 

2005 Mar. 

408.37 ±0.65 

409.398 

41.49 -41.73°S 

2 

2005 Apr. 

406.87 ± 0.76 

413.178 

2.07 - 7.30°S 

3 

2005 May 

415.06 ±6.37 

415.612 

58.86 -59.94°S 

8 

2005 Jun. 

418.52 ±4.11 

416.272 

39.31 -53.07°S 

7 

2005 Jul. 

419.72 ±0.0 

417.270 

58.00°S 

1 

2005 Aug. 

407.11 ±3.24 

417.235 

5.29- 14.23°S 

2 

2005 Oct. 

407.95 ±3.51 

413.416 

3.84 - 32.56°S 

14 

2005 Nov. 

415.00 ±4.74 

412.697 

33.3 - 58.47°S 

31 

2005 Dec. 

415.54 ±6.12 

411.300 

34.19 -56.82°S 

10 

2006 Jan. 

417.20 ±6.81 

410.560 

56.51 -59.35°S 

4 

2006 Feb. 

408.01 ±2.47 

411.206 

6.93 - 34.78°S 

17 

2006 Mar. 

414.66 ±2.85 

412.284 

37.75 -56.41°S 

8 

2006 Apr. 

410.90 ±3.32 

415.709 

9.24 - 26.55°S 

6 

2006 May 

420.28 ±7.17 

418.502 

55.85 - 
59.61°S 

9 

2006 Jun. 

421.24 ±4.56 

419.142 

49.09 - 
50.03°S 

7 

2006 Jul. 

424.89 ±5.64 

420.045 

53.13- 

59.99°S 

8 

2006 Aug. 

411.23 ±2.71 

419.937 

12.62 -25.65°S 

8 

2006 Sep. 

424.89 ± 0.0 

417.964 

59.36°S 

1 

2006 Oct. 

412.48 ±5.67 

417.042 

6.54 - 47.79°S 

20 

2006 Nov. 

417.53 ±5.90 

415.444 

29.73 - 58.47°S 

22 
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2006 Dec. 

419.54 ±6.14 

414.121 

55.19 -58.68°S 

4 

2007 Jan. 

415.44 ±4.14 

413.375 

34.8 - 58.75°S 

16 

2007 Feb. 

411.56 ±2.14 

413.885 

15.43 -32.03°S 

3 

2007 Mar. 

416.81 ±0.0 

415.333 

56.56°S 

1 

2007 May 

423.31 ±3.81 

421.244 

57.71 -59.93°S 

4 

2007 Jun. 

421.25 ±5.01 

421.572 

50.55 - 54.58°S 

4 

2007 Jul. 

429.36 ±0.0 

422.711 

52.51°S 

1 

2007 Aug. 

416.11 ±4.21 

422.476 

3.15 - 27.39°S 

5 

2007 Sep. 

421.37 ±2.91 

420.479 

50.83 -51.96°S 

3 

2007 Oct. 

414.49 ±4.20 

419.451 

0.56 - 46.54°S 

13 

2007 Nov. 

419.16 ±2.46 

418.114 

33.54 -59.55°S 

9 

2007 Dec. 

415.87 ±4.54 

416.501 

34.61 -59.69°S 

16 

2008 Jan. 

417.50 ±6.35 

416.138 

35.68 -58.72°S 

25 

2008 Feb. 

411.22 ±4.21 

416.867 

3.90 - 33.25°S 

17 

2008 Mar. 

419.04 ±3.81 

417.847 

34.95 - 57.29°S 

9 


Notes: The best-fit calculated CO? VMRs were obtained from a fit to the time series 
with Equation 1 . NMOC is the number of number of valid occultations after excluding 
cloudy and/or aerosol impacted ACE measurements and those inside or on the edge of the 
polar vortex. See text for details. 
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Figure Captions 

Figure L Fraction of remaining ACE measurements vs. latitude assuming the cloud 
rejection criterion and the microwindow at_956. 7-957.0 cm" 1 . The doted line shows the 
approximate location of the tropopause as a function of latitude based on NCEP (U.S. 
National Center for Environmental Prediction) measurements. 

Figure 2. Sample spectrum and fit obtained from analysis of CO 2 with the five selected 
microwindows. The location and adjusted tangent altitude of the measured spectrum is 
indicated at top. Below are shown the residuals on an expanded scale from the 
nonnalized measured spectrum above the measured spectrum. The tangent altitude in 
the ACE version 2.2 database of the spectrum is 6.63 1 km. 

Figure 3. Comparison of ACE monthly-mean Xco 2 at 7-10 km altitude at latitudes of 0°- 
60°N (upper panel) and 0°-60°S latitude (lower panel). No altitude shifts have been 
applied to the observations (see text for details). Open triangles and vertical lines 
indicate the monthly mean dry air mole fraction and its standard deviation, respectively. 
Measured minus calculated residuals derived with Equation 1 (solid line) are shown 
above on an expanded vertical scale. The ACE time series are compared with monthly 
average Xco 2 surface station measurements in both hemispheres. Thick lines show linear 
trend fits from the NOAA ESRL monthly mean regression fits. Colored thin lines 
connect monthly means from individual stations. Latitude, longitude, and elevation of 
NOAA ESRL stations with station codes are shown to the right. 
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Figure 4. Comparison of monthly-mean time series of ACE version 2.2 CO? monthly- 
mean Xco 2 at 7-10 km altitude between latitudes of 0° to 60°N latitude (upper panel) and 
0-60°S latitude (lower panel) with monthly-mean NOAA/ESRL surface station 
measurements (same sites station codes as in Figure 3). The ACE monthly-average 
measurements (open triangles) with a vertical line indicating the standard deviation have 
been fitted with a model that includes a sinusoidal seasonal cycle based on Equation 1 
(solid line). Altitude shifts have been applied to the measurements (see text for details). 
Measured minus calculated residuals derived with Equation 1 (solid line) are shown 
above on an expanded vertical scale. 

Figure 5 Comparisons of ACE 40°N-60°N latitude monthly average altitude-shifted 
measurements results and a best-fit to that time series with the corresponding results from 
GLOBAL VfEW model weekly time series for 7 locations near 40°N. The fonnat for the 
ACE measurements is the same as in Figures 3 and 4 and both ACE and model time series 
cover the same time period (2004 to 2008) with equation 1 used to fit the time series. The 
station code for each GLOBAL VIEW site [www.esrl.noaa.gov/gmd/ccgg/globalview/] is 
also displayed. The ACE seasonal cycle maximum occurs in the middle of January with an 
amplitude 1.158±031%. The month of seasonal cycle maximum for GLOBAL VIEW is 
found to occur in April and the amplitude is 0.92±0.01%. A total of 36 monthly means 
from ACE observations of 217 occultations are included in the analysis. 
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